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ABSTRACT 

We  present  new  observations  of  the  plerionie  supernova  remnant  3C  58  with  the  VLA  at  74  and  327 
MHz.  In  addition,  we  re-redueed  earlier  observations  at  1.4  and  4.9  GHz  taken  in  1973  and  1984.  Com¬ 
paring  these  various  images,  we  find  that  (1)  the  remnant  has  a  flat  and  relatively  uniform  speetral-index 
distribution,  (2)  any  expansion  of  the  remnant  with  time  is  significantly  less  than  that  expected  for 
uniform,  undecelerated  expansion  since  the  generally  accepted  explosion  date  in  1181  A.D.,  and  (3)  there 
is  no  evidence  for  a  nonthermal  synchrotron-emission  shell  generated  by  a  supernova  shock  wave,  with 
any  such  emission  having  a  surface  brightness  of  less  than  1  x  10“^^  W  m“^  Hz“^  sr“^  at  327  MHz. 
Subject  headings:  ISM:  individual  (3C  58)  —  radio  continuum:  ISM  —  supernova  remnants 


1.  INTRODUCTION 

Plerionie,  or  filled-center,  supernova  remnants  (Weiler  & 
Panagia  1978;  Weiler  &  Sramek  1988)  have  become  estab¬ 
lished  as  a  class.  The  supernova  remnant  3C  58 
(G130.7  -b  3.1)  was  an  early  member  of  the  class,  and  it  has 
often  been  described  as  morphologically  similar  to  the  Crab 
Nebula,  the  Crab  being  the  prototype  for  this  class.  3C  58  is 
at  a  similar  distance  as  the  Crab,  3.2  kpc  (Roberts  et  al. 
1993),  has  a  similar  size,  6'  x  9',  and  shows  a  similar  fila¬ 
mentary,  filled-center  morphology  in  the  radio  (Weiler  & 
Seielstad  1971;  Wilson  &  Weiler  1976;  Green  1986;  Rey¬ 
nolds  &  Aller  1988).  Like  the  Crab,  3C  58  has  a  relatively 
flat  radio  spectrum,  with  a  =  —  0.10  ±  0.02  (S  oc  v"^"),  from 
40  MHz  to  15  GHz  (Wilson  &  Weiler  1976;  Green  1986). 
3C  58  has  been  moderately  well  identified  with  the  histori¬ 
cal  supernova  event  of  1181  A.D.  (Clark  &  Stephenson 
1977),  which  would  give  it  an  age  similar  to  that  of  the  Crab, 
which  is  associated  with  SN  1054. 

Plerions  are  thought  to  be  powered  by  pulsars,  which 
build  their  nonthermal  nebulae  through  the  injection  of 
relativistic  particles  and  the  amplification  of  magnetic  fields, 
although  the  exact  injection  and  acceleration  processes  are 
not  yet  well  understood.  In  the  optical,  radial  velocities  less 
than  900  km  s”^  have  been  measured  for  the  filaments  in 
3C  58  (Fesen  1983),  which  is  consistent  with  the  association 
with  SN  1181  only  if  the  remnant  has  been  substantially 
decelerated  or  if  the  filaments  are  due  to  the  excitation  of  in 
situ  matter  overrun  by  the  supernova  shock.  However,  the 
linear  size  of  less  than  9  pc  (for  a  distance  of  3.2  kpc),  along 
with  the  age  of  only  820  yr,  implies  rather  low  expansion 
velocities,  which  suggests  a  low-energy  supernova  event.  In 
addition,  the  radio  flux  density  has  been  shown  to  be 
increasing  (Aller  &  Reynolds  1985;  Aller,  Aller,  &  Reynolds 
1986;  Green  1987).  It  seems  (Salvati  et  al.  1998;  Woltjer  et 
al.  1997)  that  a  nonstandard  evolution  of  the  pulsar’s  output 
is  required  to  reconcile  this  and  the  relatively  low  frequency 
of  the  observed  spectral  break  in  the  synchrotron  emission 
(Green  &  Scheuer  1992)  with  an  age  of  820  yr. 

2.  OBSERVATIONS  AND  DATA  REDUCTION 

We  observed  3C  58  using  the  A  and  B  configurations  of 
NRAO’s  Very  Large  Array  (VLA)  at  74  MHz  using  63 


spectral  channels  for  a  total  bandwidth  of  1.5  MHz,  and  at 
327  MHz  using  31  spectral  channels  for  a  total  bandwidth 
of  3  MHz.  Spectral-line  mode  was  used  to  be  able  to  excise 
radio  frequency  interference  (RFI). 

The  data  were  calibrated  and  reduced  using  NRAO’s 
AIPS  software  package.  The  instrumental  bandpass  was 
determined  using  observations  of  Cyg  A.  The  amplitude 
gains  were  set  using  Cyg  A  at  74  and  3C  48  at  327  MHz, 
and  were  tied  to  the  flux  scale  of  Baars  et  al.  (1977). 
Reduction  of  74  MHz  data  requires  somewhat  different 
techniques.^  In  particular,  RFI  is  a  significant  problem,  and 
the  data  were  first  automatically  edited  for  deviant  samples 
in  frequency  space  by  the  AIPS  task  SPFLG.  We  then  aver¬ 
aged  the  data  into  a  manageable  number  of  channels  and 
further  edited  each  channel  manually. 

Sidelobes  from  sources  well  outside  the  primary  beam  are 
a  significant  problem  at  74  MHz,  and  Cas  A,  which  was  20° 
away,  produced  strong  sidelobes  in  our  field  of  view.  In 
order  to  correctly  remove  these  sidelobes,  three- 
dimensional  imaging  was  required  (see  Perley  1999),  and  the 
AIPS  task  IMAGR  was  used  to  simultaneously  image  and 
deconvolve  numerous  flat  but  noncoplanar  regions.  These 
were  then  reassembled  into  a  single  image  by  the  task 
FLATN. 


3.  RESULTS 
3.1.  Images 

The  full-resolution  image  at  74  MHz  is  shown  in  Figure 
1.  This  image  is  made  from  the  combined  A-  and  B-array 
data  sets  with  uniform  weighting,  and  has  a  restoring  beam 
size  of  26"  (FWHM).  The  peak  brightness  for  3C  58  is  0.55 
Jy  beam”  the  background  rms  is  20  mJy  beam”  and  the 
total  flux  density  for  3C  58  is  33.6  Jy.  The  image  at  327 
MHz  is  shown  in  Figure  2.  It  is  made  from  the  combined  A- 
and  B-array  data  sets,  and  has  a  restoring  beam  size  of  8''2 
(FWHM),  a  peak  brightness  for  3C  58  of  63  mJy  beam”^,  a 
background  rms  of  0.57  mJy  beam”^,  and  a  total  flux 
density  for  3C  58  of  33.9  Jy.  It  should  be  noted  that  no 
compact  emission  source  is  detected  near  the  center  of  the 
nebula.  In  addition,  in  Figure  3,  we  show  a  low-resolution 


^  For  a  tutorial,  see  http://rsd-www.nrl.navy.mil/7213/lazio/tutorial. 
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Fig.  1. — Image  of  3C  58  at  74  MHz  from  combined  VLA  A-  and 
B-array  data.  The  FWHM  resolution  is  26"  and  is  indicated  at  lower  left, 
the  peak  brightness  is  0.55  Jy  beam"*,  and  the  background  rms  is  20  mjy 
beam"*.  The  gray  scale  is  in  mJy  beam"*.  The  contours  are  drawn  at 
-  6%,  6%,  10%,  14.1%,  20%,  28.2%,  and  40%  of  the  peak  brightness. 


image  of  the  whole  primary  beam  area  at  327  MHz,  con¬ 
volved  to  a  resolution  of  33.  No  shell  emission  is  visible. 
The  average  background  rms  between  the  edge  of  the  field 
and  the  center  is  37  mJy  beam“^,  or  3.4  x  10“^^  W  m~^ 
Hz~^  sr”^.  We  can  accordingly  place  a  3  tr  limit  on  any 
shell  emission  ofl  x  10“^^Wm“^  Hz“^  sr“^  at  327  MHz. 
While  the  field  of  view  at  74  MHz  is  even  larger,  the  sensi¬ 
tivity  is  considerably  poorer,  so  the  higher  frequency  gives 
our  most  sensitive  limit  on  extended  emission.  We  note 
however,  that  the  74  MHz  observations  would  be  sensitive 
to  very  large  shells,  i.e.,  with  diameters  greater  than  1?5,  and 
we  see  no  such  shell  emission. 

3.2.  Spectral  Index  Distribution 

We  made  an  image  of  the  spectral  index  between  74  and 
327  MHz,^  otoioT*  as  follows.  The  327  MHz  image  was  con¬ 
volved  to  the  larger  resolution  at  74  MHz,  and  then  the  two 

^  We  refer  to  speetral  indices  with  sub-  and  superscripts  indicating  the 
relevant  frequencies  in  GHz. 
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Fig.  3. — Low-resolution  image  of  the  region  surrounding  3C  58  at  327 
MHz.  The  image  has  been  convolved  with  a  Gaussian  of  33  FWHM, 
corrected  for  the  primary  beam  response,  and  blanked  where  the  primary 
beam  response  is  less  than  50%.  The  peak  flux  in  the  image  is  15.9  Jy 
beam"*.  The  background  rms  ranges  from  ~20  mjy  beam"*  near  3C  58 
to  ~45  mJy  beam"  *  near  the  edge  of  the  field  of  view.  The  gray  scale  is  in 
mJy  beam"*. 


images  were  aligned  by  the  peaks  of  the  point  sources  near 
3C  58  (at  74  MHz,  ionospheric  refraction  can  cause  sub¬ 
stantial  displacement  of  the  entire  image).  This  image  of 
aoioT,  which  has  a  resolution  of  26",  is  shown  in  Figure  4 
(regions  where  the  uncertainty  in  a®;”  is  larger  than  0.4 
have  been  blanked  for  clarity).  The  spectral  index  has  a 
mean  value  of  —0.04  over  the  nebula,  with  an  rms  of  0.18. 
Only  random  variations  are  visible  in  Figure  4.  In  particu¬ 
lar,  no  steepening  toward  the  edge  of  the  nebula  is  observed. 

We  can  calculate  the  integrated  value  of  a®;”  from  the 
total  flux  densities  in  the  naturally  weighted  images,  and  we 
obtain  a®;”  (integrated)  =  —0.07  ±  0.05,  where  the  uncer¬ 
tainty  comes  chiefly  from  the  uncertainty  in  the  VLA  ampli- 
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Fig.  2. — Image  of  3C  58  at  327  MHz  from  VLA  A-  and  B-array  data.  RIGHT  ascension  (B1950) 

The  FWHM  resolution  is  8'.'2  and  is  indicated  at  lower  left,  the  peak 

brightness  is  63  mjy  beam"*,  and  the  background  rms  is  0.57  mjy  Fig.  4. — Radio  spectral  index  of  3C  58  between  74  and  327  MHz,  ag  o), 

beam"  *.  The  gray  scale  is  in  mJy  beam"  *.  The  contours  are  drawn  at  —2,  shown  in  gray  scale.  The  FWHM  resolution  is  21"  (shown  at  lower  left). 

2, 2.8, 4, 5.7, 8, 11.3,  and  16  mJy  beam"  *.  For  reference,  we  repeat  the  total  intensity  contours  from  Fig.  1. 
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Fig.  5. — Radio  spectral  index  of  3C  58  between  327  MHz  and  4.9  GHz, 
aj  33,  shown  in  gray  scale.  The  FWHM  resolution  is  8'3  (shown  at  lower 
left).  For  reference,  the  4%,  8%,  and  16%  contours  from  the  total  intensity 
image  at  327  MHz  are  shown  also. 


tude  calibration.  This  value  is  consistent  with  ao.oa  = 

—  0.10  ±  0.02,  given  by  Green  (1986). 

We  also  reedited  and  reimaged  the  data  of  Reynolds  & 
Aller  (1988)  and  made  a  speetral-index  image  between  327 
MHz  and  4.9  GHz,  aoJs,  which  is  shown  in  Figure  5.  Again, 
no  particular  features  are  visible  in  the  body  of  the  nebula, 
and  no  steepening  toward  the  edge  is  observed.  The  mean 
and  rms  scatter  of  the  speetral  index  over  the  central  region 
of  3C  58  are  —0.051  and  0.014,  respectively.  There  is  a  very 
slight  suggestion  of  a  steepening  of  the  spectral  index  near 
the  filaments  of  less  than  0.02,  which  is  likely  not  real,  since 
systematic  sources  of  error,  in  particular  image  fidelity 
limits  (see  e.g.,  Perley  1999),  introduce  errors  in  the  spectral 
index  at  this  level  that  are  likely  eorrelated  with  the  strue- 
ture  in  the  image  itself.  The  integrated  value  of  is 

—  0.06  ±  0.03,  consistent  with  the  above  value  from  Green 
(1986). 


4.  EXPANSION  OF  3C  58 

In  order  to  determine  the  rate  of  expansion  of  the  syn¬ 
chrotron  nebula,  we  compared  our  results  with  those  from 
the  aforementioned  observations  in  1984  by  Reynolds  & 
Aller  (1988),  and  from  observations  in  1972  by  Wilson  & 
Weiler  (1976).  We  reedited  and  reimaged  both  these  sets  of 
observations  using  maximum  entropy  method  (MEM) 
deconvolution.  In  Table  1,  we  give  details  of  the  images 
used  for  determining  the  expansion.  All  the  images  at  1.4 
and  5  GHz  have  been  corrected  for  the  attenuation  of  the 


respeetive  primary  beam  patterns  (at  lower  frequencies,  the 
correetion  is  negligible).  Our  goal  is  to  determine  the  expan¬ 
sion  speed  of  the  radio  nebula.  Since  there  are  few  well- 
defined  compact  features,  we  ehoose  to  measure  the 
expansion  not  by  determining  the  proper  motion  of  individ¬ 
ual  features,  but  by  determining  an  overall  statistical  scaling 
between  our  images.  We  aceomplish  this  by  using  the 
MIRIAD  task  IMDIFF,  whieh  caleulates  how  to  make  one 
image  most  closely  resemble  another  by  calculating 
unbiased  estimators  for  the  scaling  in  size,  e,  the  scaling  and 
the  offset  in  flux  density,  A  and  b,  respectively,  and  the 
offsets  in  R.A.  and  deck,  x  and  y,  respectively,  needed  to 
make  the  second  image  most  closely  resemble  the  first.  We 
are  chiefiy  interested  in  the  expansion  factor  e,  but  because 
of  uneertainties  in  flux  calibration,  absolute  position,  and 
image  zero-point  offsets  caused  by  missing  short  spacings, 
we  determine  all  parameters.  This  method  was  developed 
by  Tan  &  Gull  (1985),  and  more  details  of  its  use  in  a  similar 
situation  are  given  by  Bietenholz  et  al.  (1991). 

Since  we  have  shown  above  that  the  radio  spectral  index 
is  uniform  to  within  the  uneertainties,  images  at  any  two 
frequeneies  and  V2  should  closely  resemble  one  another, 
except  for  a  flux-scaling  factor,  which  is  given  by  (V1/V2)'', 
where  a  is  the  spectral  index.  Thus,  we  ean  usefully  compare 
images  taken  at  different  observing  frequencies  to  test  for 
expansion. 

We  ran  IMDIFF  on  various  combinations  of  the  above 
images,  always  with  both  images  convolved  to  the  same 
effeetive  resolution,  to  get  the  expansion  factors  shown  in 
Table  2.  In  order  to  estimate  the  uncertainty  in  our  calcu¬ 
lation  of  e,  we  also  ran  IMDIFF  on  two  pairs  of  (almost) 
contemporaneous  images,  for  which  e  should  be  0.0.  They 
are  the  expansion  between  1.4  and  4.9  GHz  at  essentially 
the  same  epoch,  and  the  expansion  between  a  CLEAN  and 
a  MEM  image  of  the  1973  4.9  GHz  data,  and  the  results  are 
shown  in  the  last  two  lines  in  Table  2.  In  both  of  these  test 
cases,  we  obtain  expansions  of  less  than  0.2%,  which  sug¬ 
gests  that  the  uncertainties  in  our  expansion  factors  are  of 
this  order. 

In  order  to  determine  the  expansion  speed,  we  fit  a 
straight  line  through  the  first  seven  expansion  factors  in 
Table  2  by  weighted  least  squares.  The  expansion  for  zero 
time  difference  must  be  0.0,  so  we  constrain  the  fit  line  to 
pass  through  the  origin.  Sinee  we  do  not  have  formal  uncer¬ 
tainties  in  the  expansion  factor,  we  use  the  weights  listed  in 
Table  2,  which  we  derived  as  follows:  it  seems  reasonable  to 
suppose  that  the  uncertainty  in  the  expansion  factor  is  pro¬ 
portional  to  the  rms  of  the  residuals  of  the  IMDIFF  fit,  (7^, 
divided  by  the  average  brightness  in  the  images — more  fiux 


TABLE  1 

Images  Used  for  the  Expansion  Determination 


Frequency  Image  rms  Convolving  Beam 

Date  Telescope  (MHz)  (Jy)  (mjy  per  beam)  (arcsec  x  arcsec,  at  deg)  Reference 


1973  Feb-Jun  .  WSRT  4995  30.1  2.44  8.3  x  8.3  1 

1973  Dec  15 .  WSRT  1415  33.1  1.35  22.7  x  20.9,  at  0.0  1 

1984  Jan-Aug .  VLA  4866  29.5  0.43  8.3  x  8.3  2 

VLA  4866  29.5  1.71  17.7  x  17.4,  at  -7.5  2 

1984  Jan-Dec  .  VLA  1446  33.9  0.49  17.7  x  17.4,  at  -7.5  2 

1998  Oct  1 .  VLA  327  33.0  1.17  17.7  x  17.4,  at  -7.5  3 

1998  Mar-Oct .  VLA  327  34.0  0.58  8.3  x  8.3  4 


“  is  the  total  cleaned  flux  density  in  the  image. 

References. — (1)  Wilson  &  Weiler  1976;  (2)  Reynolds  &  Aller  1988;  (3)  This  paper,  B-array  only;  (4)  This  paper,  A  and  B 
arrays 
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TABLE  2 

Expansion  of  3C  58 


First  Image  Second  Image 


Year* 

Frequency 

(MHz) 

Year* 

Frequency 

(MHz) 

Resolution'’ 

(arcsec) 

Time  Interval 
(yr) 

Change  in  Size 
(%) 

rms  of  Fit 
(mjy  beam"') 

Weight' 

1984.39 

1.4 

1998.75 

0.3 

17 

14.4 

0.80 

1.4 

1.1 

1984.27 

4.9 

1998.46 

0.3 

8 

14.2 

0.45 

0.6 

1.7 

1973.96 

1.4 

1998.75 

0.3 

21 

24.8 

0.80 

3.0 

0.5 

1973.96 

1.4 

1984.39 

1.4 

21 

10.4 

-0.10 

2.2 

0.9 

1973.96“ 

1.4 

1984.39“ 

1.4 

21 

10.4 

0.05 

1.8 

1.3 

1973.29“ 

4.9 

1998.46“ 

0.3 

8 

25.2 

0.15 

0.9 

0.8 

1973.29“ 

4.9 

1984.27“ 

4.9 

8 

11.0 

-0.10 

0.7 

1.3 

The  results  below  are  included  as  a  check 

1984.27 

4.9 

1984.39 

1.4 

17 

0.1 

-0.20 

1.4 

1.4 

1973.27“ 

4.9 

1973.27' 

4.9 

8 

0.0 
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“For  observations  involving  multiple  arrays,  the  date  given  is  the  average  over  all  the  arrays,  weighted  by  the  length  of  the  individual  array 
observing  runs. 

’’  The  rounded  average  resolution,  see  Table  1. 

“  These  are  the  weights  used  for  fitting  the  overall  expansion;  see  text  for  details. 

‘‘  These  images  have  been  high-pass  filtered  using  a  circular  Gaussian  of  FWHM  60". 

'  This  line  pertains  to  the  difference  between  CLEAN  and  MEM  images  from  the  same  data 


in  the  image  allows  a  relatively  better  fit  for  the  same  rms. 
Since  the  total  flux  density  is  roughly  the  same  in  all  the 
images,  the  average  brightness  is  inversely  proportional  to 
the  number  of  beam  areas  over  the  source.  It  also  seems 
reasonable  that  the  uncertainty  is  inversely  proportional  to 
the  square  root  of  the  number  of  beam  areas  in  the  fitting 
region,  We  thus  adopt  as  the  weight 


W  = 


1 


4,  A 


fitting  region 


(1) 


For  the  5  GHz,  1973  Westerbork  Synthesis  Radio  Tele¬ 
scope  (WSRT)  data,  we  use  only  expansions  derived  from  a 
version  of  the  image  high-pass  filtered  with  a  Gaussian  of  1' 
FWHM,  since  this  image  has  evident  short-spacing  prob¬ 
lems  (Wilson  &  Weiler  1976).  The  points,  as  well  as  the 
weighted  least-squares  fit,  are  plotted  in  Figure  6. 

The  straight-line  fit  to  the  data  in  Table  2  by  weighted 
least  squares  is  plotted  as  the  solid  line  in  Figure  6.  The 
slope,  or  average  expansion  speed,  that  we  obtain  from  the 
fit  is  (0.020  +  0.008)%  yr“^,  which  corresponds  to 
-^910  ±  360  and  550  ±  220  km  s“^  along  the  major  and 
minor  axes  of  the  nebula,  respectively.  We  can  estimate  the 
uncertainties  from  the  scatter  of  the  data  points  about  the  fit 
line,  and  we  find  that  a  W  =  I  corresponds  to  an  uncer¬ 
tainty  of  0.34%  in  the  expansion  factor.  Our  two  check 
points,  which  measure  the  expansion  between  almost  simul¬ 
taneous  images  (the  last  two  lines  of  Table  2),  are  indeed 
consistent  with  an  expansion  of  0.0,  as  they  should  be.  We 
also  note  that  an  unweighted  fit  gives  almost  the  same 
expansion  speed  and  uncertainty. 

Our  measured  expansion  speed  (solid  line  in  Fig.  6)  is 
much  less  than  would  be  expected  for  an  undecelerated 
expansion  of  3C  58  since  1181  A.D.  The  undecelerated 
“historical”  expansion  speed  would  be  0.124%  yr“^  and  is 
plotted  as  the  dashed  line  in  Figure  6.  In  addition,  in  Figure 
7,  we  show  profiles  drawn  through  the  long  axis  of  3C  58  at 
1.4  GHz  in  1984  and  at  0.3  GHz  in  1998,  and  for  compari¬ 
son,  a  hypothetical  profile  that  is  generated  by  unde¬ 
celerated  expansion  of  the  1984  profile  to  1998.  This  last  can 


be  seen  to  be  notably  larger  than  the  true  1998  profile,  also 
indicating  a  lower  than  “  historical  ”  expansion  speed. 

The  nebular  structure  is  dominated  by  the  relatively 
sharp  outer  edges  rather  than  by  any  well-defined  filaments 
or  clumps,  so  our  expansion  results  derive  in  considerable 
part  from  the  overall  extent  of  the  nebula.  Unfortunately, 
this  makes  our  expansion  results  sensitive  to  short-spacing 
problems  and  also  to  the  accuracy  of  the  beam  corrections, 
especially  at  5  GHz,  where  the  VLA  and  WSRT  primary 
beams  are  825"  and  650"  FWHM,  respectively.  Because  of 
this,  we  cannot  exclude  the  presence  of  systematic  biases  in 
our  results.  However,  we  feel  that  the  aggregate  expansion 
speed,  which  is  derived  from  a  variety  of  combinations  of 


Fig.  6. — Expansion  of  3C  58  over  various  time  intervals.  The  points 
plotted  represent  the  expansion  required  by  IMDIFF  over  the  pairs  of 
observations  from  Table  2,  excluding  the  last  two  check  points.  The  soUd 
line  indicates  the  fit  expansion  speed  of  (0.020  +  0.008)%  yr“\  and  the 
dashed  line  indicates  the  expansion  of  0.124%  yr“^  expected  for  unde¬ 
celerated  expansion  since  1181  A.D. 
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Fig.  7. — Profiles  through  3C  58.  The  top  panel  shows  the  profile 
through  3C  58,  drawn  at  P.A.  —86°,  i.e.,  through  the  long  axis  of  the 
nebula.  The  solid  line  represents  the  flux  density  at  1.4  GHz  in  1984.  The 
dashed  line  represents  the  flux  density  at  0.3  GHz  in  1998,  multipUed  by  a 
flux  density  scaling  factor  so  as  to  most  closely  match  the  1984  profile  (the 
scaling  factor  was  determined  by  IMDIFF,  see  text  §  4).  The  dotted  Une, 
drawn  for  comparison,  represents  the  1984  data  homologously  expanded 
by  1.74%,  which  is  the  expected  amount  if  3C  58  is  the  undecelerated 
remnant  of  SN  1181.  The  middle  and  the  bottom  panel  are  magnified 
regions  of  the  profile  in  the  top  panel,  showing  details  of  the  eastern  and 
western  limbs  of  3C  58,  respectively.  The  FWHM  resolution  is  17'.'7 
X  17'.'4,at  -7?5. 

different  epochs,  resolutions,  telescopes,  and  frequencies,  is 
unlikely  to  be  greatly  biased. 

5.  DISCUSSION 

We  measure  an  expansion  rate  of  the  radio  nebula  for 
3C  58  of  (0.020  +  0.008)%  yr“^  (corresponding  to 
-^910  ±  360  and  550  ±  220  km  s“^  along  the  major  and 
minor  axes  of  the  nebula,  respectively).  This  is  much  smaller 
than  would  be  expected  if  3C  58  were  indeed  the  unde¬ 
celerated  remnant  of  SN  1181,  which  would  require  an 
expansion  of  0.124%  yr“^  (~5500  and  ^3300  km  s“^ 
along  the  major  and  minor  axes,  respectively).  Conversely, 
our  measured  expansion  rate  suggests  an  age  for  3C  58  of 
5000  yr  with  a  3  a  lower  limit  of  2250  yr,  which  is  incompat¬ 
ible  with  the  generally  accepted  association  with  SN  1181. 
Although  we  cannot  entirely  rule  out  systematic  sources  of 
error  in  our  determination  of  the  expansion,  we  think  it 
unlikely  that  they  are  large  enough  to  make  our  results 
compatible  with  an  expansion  age  of  ~  820  yr. 

This  low  expansion  speed  is  corroborated  by  the  proper 
motions  of  optical  filaments:  Fesen,  Kirshner,  &  Becker 
(1988)  and  van  den  Bergh  (1990)  looked  for  such  proper 
motions  and  failed  to  find  any  larger  than  0'(07  yr“^.  Fesen 


(1983)  measured  optical  radial  velocities  for  several  fila¬ 
ments  in  3C  58.  The  largest  radial  velocities  he  measured 
were  ~900  km  s“^  for  filaments  situated  near  the  projected 
center  of  the  nebula,  where  the  expansion  is  along  the  line  of 
sight,  and  thus  where  the  expansion  speed  is  expected  to  be 
approximately  equal  to  the  radial  velocity.  This  is  also  sub¬ 
stantially  lower  than  the  >3000  km  s“^  implied  by  unde¬ 
celerated  expansion  and  a  distance  of  3.2  kpc.  (We  note  that 
Fesen,  who  used  a  somewhat  smaller  distance  to  3C  58  of 
2.6  kpc,  reconciled  these  radial  velocities  with  an  age  of  820 
yr  by  suggesting  that  substantial  deceleration  had  taken 
place,  and  that  the  average  speed  since  1181  A.D.  was 
~3000  km  s“^.)  In  order  to  determine  an  expansion  age 
from  this  radial  velocity,  we  need  to  know  the  radius  at 
which  the  measured  filament  lies  from  the  center  of  the 
nebula.  If  we  take  this  radius  to  be  the  mean  projected 
radius  of  ~4  pc,  these  radial  velocities  also  imply  a  large 
age  of  -^4400  yr,  but  since  the  optical  emission  is  faint, 
those  radial  velocities  may  have  been  measured  for  fila¬ 
ments  nearer  the  center  of  the  nebula,  and  thus  a  smaller 
radius  and  consequently  a  smaller  age  is  possible.  However, 
we  consider  it  unlikely  that  the  largest  radial  velocities  were 
measured  for  filaments  only  1  pc  from  the  center  of  the 
nebula.  If  3C  58  is  associated  with  SN  1181,  all  this  evidence 
suggests  that  it  must  have  been  greatly  decelerated. 

The  simplest  solution  to  these  inconsistencies  is  to 
suggest  that  3C  58  is  not  the  remnant  of  SN  1181.  This  also 
alleviates  the  difficulties  in  accounting  for  the  very  sharp 
and  rather  low  frequency  spectral  break  observed  in  3C  58 
(Woltjer  et  al.  1997;  Green  &  Scheuer  1992).  Note  that  a 
lower  value  for  the  somewhat  uncertain  distance  estimate 
for  3C  58  would  reduce  the  age  determined  from  the  radial 
velocities,  but  not  that  determined  from  the  proper  motions, 
so  that  is  not  a  solution.  The  positional  coincidence 
between  SN  1181  and  3C  58,  however,  was  confirmed  by  a 
recent  reexamination  (Stephenson  &  Green  1999),  and  so  if 
3C  58  is  not  to  be  associated  with  SN  1181,  one  must  ask 
where  the  remnant  of  SN  1181  is,  since  there  are  no  other 
obvious  candidates,  and  since  it  seems  unlikely  that  the 
remnant  would  already  be  too  faint  to  observe  at  an  age  of 
only  820  yr. 

Van  den  Bergh  (1990)  suggests  that  the  presently  visible 
optical  filaments  predate  the  supernova  event,  and  are 
quasi-stationary  but  energized  by  the  passing  supernova 
shock  wave.  It  is  not  clear  that  the  low  expansion  rate  of  the 
synchrotron  bubble  can  be  accounted  for  in  this  way,  since 
it  seems  unlikely  that  these  preexisting  filaments,  which 
have  a  low  filling  factor,  would  be  able  to  confine  the  syn¬ 
chrotron  bubble.  While  some  coupling  between  the  syn¬ 
chrotron  bubble  and  the  thermal  filaments  would  be 
expected,  since  the  magnetic  field  in  the  bubble  would  wrap 
around  the  thermal  filaments,  one  would  expect  the  syn¬ 
chrotron  bubble  to  rapidly  expand  beyond  the  filaments, 
along  with  the  rest  of  the  supernova  ejecta.  The  outside 
edge  of  the  visible  synchrotron  bubble  should  be  expanding 
at  least  at  the  speed  of  the  local  ejecta,  which  would  be 
~4000  km  s”^  for  undecelerated  expansion.  Should  preex¬ 
isting  filaments  nonetheless  be  somehow  confining  the  syn¬ 
chrotron  bubble,  then  the  filaments  should  be  accelerated. 
The  current  measurements,  unfortunately,  do  not  have  the 
precision  necessary  to  determine  if  this  is  the  case. 

A  possible  explanation  for  slow-moving  filaments  is  the 
following:  the  magnetic  field  wraps  around  the  preexisting 
filaments  and  causes  an  intensification  of  the  synchrotron 


No.  2,  2001 


RADIO  SPECTRAL  INDEX  AND  EXPANSION  OF  3C  58 


111 


emissivity  in  their  vicinity,  which  is  why  these  filaments  are 
visible  in  radio  synchrotron  emission.  If  the  structure  of  the 
radio  nebula  is  dominated  by  these  irregularities,  then  our 
determination  of  the  expansion  would  determine  the  expan¬ 
sion  precisely  of  the  filaments  and  not  of  the  synchrotron 
bubble  as  a  whole.  We  note  that  expansion  speeds  that  we 
measure  at  lower  elfective  resolutions  tend  to  be  slightly 
smaller  than  those  at  larger  resolutions,  as  would  be 
expected  in  this  case,  but  even  the  measured  expansions  at 
lower  resolutions  are  far  below  the  expected  values  for 
undecelerated  expansion.  Furthermore,  profiles  through  the 
images  do  not  support  the  conjecture  of  quasi-stationary 
filaments  in  an  expanding  bubble :  the  profiles  displayed  in 
Figure  7  show  that  even  the  expansion  of  the  edge  of  the 
nebula  appears  to  be  considerably  less  than  is  expected 
from  undecelerated  expansion. 

Let  us  turn  to  some  general  physical  considerations 
about  the  synchrotron  bubble.  The  pressure  in  the  synchro¬ 
tron  bubble  is  relatively  high:  standard  arguments  give  a 
minimum  pressure  of  greater  than  10“^°  dyn  cm“^.  This  is 
the  minimum  pressure  of  the  synchrotron  fluid  assuming  a 
filling  factor  of  1.0  and  a  line-of-sight  depth  equal  to  the 
long  axis  of  the  projected  nebula.  Both  a  lower  filling  factor 
and  the  admixture  of  any  thermal  material  would  raise  the 
pressure,  and  we  think  it  unlikely  that  the  line-of-sight 
dimension  is  much  larger  than  the  long  axis  of  the  projected 
nebula. 

The  slow  expansion  of  the  bubble  therefore  necessitates 
some  confining  pressure,  since  the  sound  and  Alfven  veloci¬ 
ties  in  the  bubble  are  also  high.  This  leads  to  the  question  of 
what  surrounds  the  presently  visible  nebula.  The  three  pos¬ 
sibilities  are  (1)  the  ISM,  (2)  the  wind  bubble  of  the  presu¬ 
pernova  progenitor  of  SN  1181,  or  (3)  a  freely  expanding 
stellar  envelope  thrown  off  in  the  supernova  explosion.  This 
last  might  be  expected,  since  the  total  energy  in  the  visible 
remnant  is  far  below  the  expected  energy  of  10®^  ergs  rel¬ 
eased  in  a  supernova  for  any  reasonable  mass  in  the  fila¬ 
ments.  This  last  situation  does  seem  to  obtain  in  the  Crab 
Nebula,  where  strong  circumstantial  evidence  points  to  the 
existence  of  such  a  freely  expanding  stellar  envelope  (Hester 
et  al.  1996;  Sankrit  &  Hester  1997),  even  though  the  radio 
emission  from  this  shell  has  so  far  eluded  detection 
(Bietenholz  et  al.  1997;  Frail  et  al.  1995),  and  the  direct 
optical  detection  is  controversial. 

The  static  pressure  of  any  of  these  three  exterior  media, 
is  expected  to  be  considerably  smaller  than  the 
minimum  pressure  in  the  synchrotron  bubble; 

Pext  =  1-4  X  10“^^  dyn  cm“^(n/l  cm“^)(r/10"^  K) ,  (2) 

where  n  is  the  number  density  and  T  is  the  temperature. 
This  suggests  that  the  synchrotron  bubble  is  confined  by 
ram  pressure  as  it  expands  into  the  surrounding  material. 
We  are  then  brought  to  a  crucial  difference  between  3C  58 
and  the  Crab  Nebula:  in  the  case  of  the  Crab,  the  synchro¬ 
tron  bubble  has  been  accelerated  since  the  supernova  event 
(Bietenholz  et  al.  1991),  and  is  therefore  expanding  into  the 
stellar  envelope,  and  confined  by  the  ram  pressure  of  this 
expansion.  In  the  case  of  3 C  58,  by  contrast,  the  synchro¬ 
tron  bubble  is  expanding  considerably  more  slowly  than 
any  freely  expanding  material  thrown  off  in  SN  1181,  and 
cannot  therefore  be  confined  by  ram  pressure  in  any  freely 
expanding  ejecta. 

This  implies  for  3C  58  that  there  is  no  freely  expanding 
envelope,  making  it  a  very  low-energy  supernova  event,  or 


that  the  envelope  has  been  greatly  decelerated,  which  seems 
unlikely  in  light  of  the  nondetection  of  any  shell  emission 
exterior  to  the  synchrotron  remnant,  or  that  it  is  not  the 
remnant  of  SN  1181.  The  ram  pressure  required  suggests 
that  the  external  medium  have  a  number  density,  n^^f,  of 

«ext  ~  0.1  cm (i;,h„ek/1000  km  s"i)^(m/mH) ,  (3) 

where  fshock  is  the  difference  in  velocity  between  the  syn¬ 
chrotron  bubble  and  the  undisturbed  exterior  medium,  and 
m  is  the  mean  mass  per  particle. 

A  possible  scenario  would  be  that  3C  58  is  the  composite 
remnant  of  a  binary  where  both  stars  have  undergone 
supernova  explosions,  with  the  earlier  supernova  leaving 
the  slowly  expanding  remnant,  and  the  later,  low-energy 
one  providing  the  historical  supernova  event  in  1181  A.D. 
This  scenario  seems  rather  contrived,  since  having  two 
supernovae  within  a  few  thousand  years  of  one  another  is 
quite  unlikely. 

We  turn  to  the  nondetection  of  any  shell  emission.  If 
there  were  an  envelope  ejected  in  SN  1181,  as  discussed 
above,  it  would  have  been  decelerated  by  at  least  a  factor  of 
3,  which  in  turn  suggests  that  it  would  have  already  swept 
up  considerable  external  material,  and  be  well  into  the 
Sedov  phase  of  its  evolution.  It  is  then  surprising  that  the 
radio  brightness  of  this  exterior  shell  is  considerably  less 
than  that  of  the  faintest  known  shells.  Again,  it  is  instructive 
to  point  out  the  difference  between  3C  58  and  the  Crab 
Nebula.  Also,  for  the  Crab,  radio  emission  from  the  shell 
has  not  yet  been  detected  to  similar  limits  in  surface  bright¬ 
ness,  but  the  envelope  is  still  in  the  free-expansion  phase 
and  is  not  expected  to  have  converted  much  of  its  energy 
into  relativistic  particles,  and  therefore  is  expected  to  be 
fainter  than  most  shell  remnants. 

Could  a  shell  for  3C  58  be  masked  because  it  is  largely 
coincident  with  the  synchrotron  nebula?  Shock-accelerated 
radio-synchrotron  emitting  electrons  typically  produce 
radio  emission  with  a  spectral  index  a  ^  —0.6.  Our  spectral 
index  results,  however,  indicate  that  there  is  no  significant 
radio  emission  with  a  =  —0.6  near  the  edge  of  the  nebula, 
since  we  observe  no  steepening  of  the  spectrum  there.  Also, 
the  total  flux  density  that  we  measure  at  74  MHz  gives  an 
integrated  spectral  index  no  steeper  than  that  at  higher  fre¬ 
quencies,  contrary  to  what  would  be  expected  if  there  were  a 
steep  spectrum  component  present. 

Both  ao;o7  and  aoJs  show  no  spatial  variation  over  the 
nebula,  showing  rms  variations  of  only  less  than  0.05. 
The  uniformity  of  the  spectral  index  suggests  that  a  single 
acceleration  mechanism  is  responsible  for  all  the  radio- 
emitting  electrons.  This  is  thought  to  be  an  as  yet  unde¬ 
tected  pulsar,  although  there  are  a  number  of  problems  with 
3C  58  in  this  regard.  Recent  X-ray  results  by  Bocchino  et  al. 
(2001)  do  not  show  the  expected  blackbody  X-ray  emission 
from  the  neutron  star.  No  pulsed  emission  has  so  far  been 
found  at  any  wavelength,  although  in  a  deep-radio  image. 
Frail  &  Moffett  (1993)  detected  a  “  wisplike  ”  feature,  similar 
to  ones  usually  associated  with  the  termination  shock  of  a 
pulsar  wind. 

6.  A  SHELL  AROUND  3C  58 

In  the  case  of  the  Crab  Nebula,  there  is  substantially  less 
kinetic  energy  in  the  remnant  than  the  10^  ^  ergs  expected 
from  a  Type  II  supernova.  This  has  suggested  the  existence 
of  a  rapidly  moving  shell  around  the  currently  visible 
remnant  that  carries  the  excess  energy  and  also  confines  the 
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relatively  high-pressure  synehrotron  bubble.  Much  obser¬ 
vational  effort  has  been  expended  in  search  of  this  shell 
around  the  Crab,  but  no  shell  has  yet  been  detected  in  the 
radio,  ruling  out  any  radio  emission  with  a  surface  bright¬ 
ness  much  less  than  that  of  the  faintest  known  young  shell 
remnant,  SN  1006  (Frail  et  al.  1995;  Bietenholz  et  al.  1997). 
However,  in  the  case  of  the  Crab,  circumstantial  evidence 
for  the  existence  of  expanding  material  outside  the  currently 
visible  nebula  has  been  made  from  Hubble  Space  Telescope 
observations  (Hester  et  al.  1996),  which  show  the  Rayleigh- 
Taylor  unstable  interface  between  the  high-pressure  syn¬ 
chrotron  bubble  and  the  hydrogen  shell. 

In  the  case  of  3C  58,  there  is  a  similar  deficit  in  energy, 
especially  in  view  of  the  low  expansion  velocity  of  both  the 
synchrotron  bubble  and  the  optical  filaments,  and  so  a 
rapidly  moving  shell  might  also  be  expected  exterior  to  the 
currently  visible  nebula.  Several  searches  for  a  radio  shell 
around  3C  58  have  been  made,  notably  by  Reynolds  & 
Aller  (1985),  and  no  sign  of  radio  emission  from  a  shell  has 
yet  been  detected. 

Our  327  MHz  results  give  a  sensitive  upper  limit  to  the 
radio  emission  from  a  putative  shell  oflxlO“^^Wm“^ 
Hz“^  sr“^  at  that  frequency.  Scaling  with  an  assumed  shell 
spectral  index  of  —0.6,  this  is  equivalent  to  4  x  10“^^  W 
m”^  Hz“^  sr“^  at  1.4  GHz,  thus  very  slightly  more  sensi¬ 
tive  than  the  earlier  limit  of  4J  x  10“^^  W  m“^  Hz“^  sr“^ 
obtained  by  Reynolds  &  Aller  (1985).  We  are,  however, 
more  sensitive  to  steep-spectrum  shells  and  also  to  larger 
shells  because  of  our  larger  field  of  view  at  327  MHz.  We 
note  that,  even  on  our  74  MHz  images,  where  the  primary 
beam  radius  is  ^  5°,  there  is  no  sign  of  shell  emission,  albeit 
at  a  lower  sensitivity  than  at  327  MHz.  A  sensitivity  to 
larger  shells  may  be  germane  if  3C  58  is,  in  fact,  older  than 
SN  1181,  since  at  10,000  km  s“^,  a  3500  yr  old  shell  at  3.2 
kpc  would  have  an  angular  diameter  of  1?2. 

7.  CONCLUSIONS 

New  observations  at  74  and  327  MHz  and  reanalysis  of 
older  observations  dating  back  to  1974  indicate  that: 

1.  3C  58  has  a  flat  spectral  index  of  =  —0.05,  with 
little  variation  over  the  nebula. 

2.  The  expansion  of  the  synchrotron  nebula  is  much 
slower  than  expected  for  an  explosion  date  in  1181  A.D., 
with  constant  undecelerated  expansion  since  then.  A  pos¬ 
sible  explanation  is  that  3C  58  is  not  the  remnant  of  SN 
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1181  but  is  far  older,  with  an  age  of  5000  yr;  however,  we 
are  reluctant  to  abandon  the  generally  accepted  association 
until  a  better  understanding  of  the  dynamics  of  plerions  is 
available,  and  a  better  candidate  remnant  for  SN  1181  is 
found.  An  alternate  explanation  is  that  3C  58  has  been 
greatly  decelerated. 

3.  As  in  the  Crab  Nebula,  any  nonthermal  emission  from 
a  shock  wave  outside  the  observable  plerion  has  a  very  low 
surface  brightness  of  less  than  1  x  10“^^  W  m~^  Hz“^ 
sr“^.  If  3C  58  has  been  greatly  decelerated  as  suggested  in 
item  2  above,  then  one  would  expect  radio  emission  from 
such  a  shock  wave. 

We  close  by  summarizing  several  apparent  inconsis¬ 
tencies  in  the  observational  material  on  3C  58.  First,  the 
measured  low  proper  motions,  both  in  the  radio  and  in  the 
optical,  along  with  the  low  synchrotron  break  frequency, 
suggest  an  age  much  older  than  820  yr.  On  the  other  hand, 
the  positional  association  with  SN  1181  seems  secure,  and 
no  other  remnant  has  been  identified  for  this  historical 
supernova.  The  radio  spectral  index  is  very  uniform,  sug¬ 
gesting  a  common  source  for  the  radio-emitting  electrons, 
namely  a  pulsar  wind,  yet  the  postulated  pulsar  has  escaped 
detection  and  the  expected  X-ray  emission  from  the  hot 
surface  of  the  neutron  star  has  not  been  detected.  Finally, 
the  currently  observable  remnant  contains  only  a  small 
fraction  of  the  kinetic  energy  expected  from  a  Type  II  super¬ 
nova.  Unless  SN  1181  was  an  anomalously  low  energy 
event,  a  shell  consisting  of  the  rapidly  moving  external 
layers  of  the  star  would  be  expected  to  carry  the  balance  of 
the  kinetic  energy.  Such  a  shell  has  not  been  detected,  and 
would  be  anomalously  faint  in  comparison  to  other  known 
shells.  We  hope  that  future  study  will  allow  us  a  more  com¬ 
plete  and  consistent  picture  of  this  fascinating  object. 
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